Self-avoiding walks and connective constants in small-world networks 
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Long-distance characteristics of small-world networks have been studied by means of self-avoiding 
walks (SAW's). We consider networks generated by rewiring links in one- and two-dimensional 
regular lattices. The number of SAW's u n was obtained from numerical simulations as a function 
of the number of steps n on the considered networks. The so-called connective constant, /j, — 
lim„_ 00 Un/un-i, which characterizes the long-distance behavior of the walks, increases continuously 
with disorder strength (or rewiring probability, p). For small p, one has a linear relation /i = jio + ap, 
Ho and a being constants dependent on the underlying lattice. Close to p = 1 one finds the behavior 
expected for random graphs. An analytical approach is given to account for the results derived from 
numerical simulations. Both methods yield results agreeing with each other for small p, and differ 
for p close to 1, because of the different connectivity distributions resulting in both cases. 
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I. INTRODUCTION 

Our world is formed by networks of different types (so- 
cial, biological, technological, economic), whose charac- 
terization has launched in last years the emergence of 
models incorporating the basic ingredients of real- life net- 
works [1-3]. In particular, social networks form the sub- 
strate where processes such as information spreading or 
disease propagation take place. One expects that the 
structure of these complex networks will play an impor- 
tant role in such dynamical processes, which are usu- 
ally studied by means of stochastic dynamics and ran- 
dom walks. Some processes, such as navigation and ex- 
ploratory behavior are neither purely random nor totally 
deterministic, and can be also described by walks on 
graphs [4]. In this context, the generic properties of de- 
terministic navigation [5] and directed self-avoiding walks 
[6] in random networks have been analyzed recently. 

In last years, networks displaying the "small- world" 
effect have been intensively studied [7-11]. Watts and 
Strogatz [7,12] proposed for this kind of networks a model 
based on a locally connected regular lattice, in which a 
fraction p of the links between nearest-neighbor sites are 
replaced by new random connections, thus creating long- 
range "shortcuts" . Hence, one has in the same network 
a local neighborhood (as for regular lattices) and some 
global properties of random graphs [13]. The small-world 
effect is usually measured by the scaling behavior of the 
characteristic path length £, defined as the average of the 
distance between any two sites. In small- world networks, 
one has a logarithmic increase of £ with the network size, 
as happens for random graphs [2,13,14]. 

This short global length scale changes strongly the 
behavior of statistical physical problems on small-world 
networks, as compared with regular lattices (where one 
has I ~ N 1 ^, N being the system size and d the lat- 



tice dimension). Among these problems, one finds signal 
propagation [7], spread of infections [15,16], and random 
spreading of information [17-19]. Site and bond perco- 
lation [16,20], as well as the Ising [21-23] and XY mod- 
els [24], have been also studied in these networks. Most 
of the published work on small worlds has focussed on 
networks obtained from one-dimensional lattices (rings). 
Small- world networks built by rewiring lattices of higher 
dimensions have being employed to study percolation, as 
a model of disease propagation [20,25] . Several character- 
istics of random walks on this kind of networks have been 
analyzed in connection with diffusion processes [26,27]. 
In particular, some properties of these walks, such as the 
probability of returning to the origin, were found to be 
intermediate between those corresponding to fractals and 
Cay ley trees [28]. 

In this paper we study self-avoiding walks in small- 
world networks built up from one- and two-dimensional 
regular lattices. A self-avoiding walk (SAW) is defined 
as a walk along the bonds of a given network which can 
never intersect itself. The walk is restricted to moving 
to a nearest-neighbor site during each step, and the self- 
avoiding condition constrains the walk to occupy only 
sites which have not been previously visited in the same 
walk [29] . SAW's have been used for modeling the large- 
scale properties of long-flexible macromolecules in solu- 
tion [30], as well as for the study of polymers trapped 
in porous media, gel electrophoresis, and size exclusion 
chromatography, which deal with the transport of poly- 
mers through membranes with small pores [31]. They 
have been also employed to characterize complex crystal 
structures [32] and to analyze critical phenomena in lat- 
tice models [29,33]. Universal constants for SAW's have 
been discussed by Privman et al. [34]. 

The paper is organized as follows. In Sec. II we give 
some basic definitions and concepts related to SAW's. In 
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Sec. Ill we present results for SAW's on simulated small- 
world networks, and in Sec. IV we give an approximate 
analytical procedure to calculate the number of SAW's 
on this kind of networks. The paper closes with some 
conclusions in Sec. V. 



II. BASIC DEFINITIONS 

For regular lattices, the number u n of different SAW's 
starting from a generic site has an asymptotic depen- 
dence for large n: [34,35] 



TTl 

Qrd(m) = — Prd(m) 

z 



(4) 



u„ ~ n 



(1) 



where 7 is a critical exponent which depends on the lat- 
tice dimension, and \x is the so-called "connective con- 
stant" or effective coordination number of the considered 
lattice [35,36]. In general, for a lattice with coordination 
number (or connectivity) z, one has zz < z — 1. This 
parameter /x can be obtained from Eq. (1) by the limit: 



v. 



lim 



(2) 



The connective constant depends upon the particu- 
lar topology of each lattice, and has been determined 
very accurately for two-dimensional (2D) and three- 
dimensional (3D) lattices [37]. In the following, we will 
consider Eq. (2) as a definition of the connective constant 
fi for any network. [Note that the limit in Eq. (2) is well 
defined provided that the mean connectivity is finite.] 

For random and small-world networks the number of 
SAW's of length n depends on the considered starting 
node of the network. In the sequel we will call u n the 
average number of SAW's of length n, i.e. the mean 
value obtained (for each n) by averaging over the net- 
work sites. For small- world networks one expects \x val- 
ues larger than that corresponding to the starting regu- 
lar lattice. In particular, \x is expected to increase with p 
and approach the value corresponding to random lattices 
with mean connectivity z as p — » 1. 

For regular lattices all nodes have the same connectiv- 
ity, i.e., the same number of nearest neighbors. However, 
for p > different connectivities m are possible, giving 
rise to a probability distribution for which analytic ex- 
pressions have been found [19,21]. For a (large) random 
network with mean connectivity z, the connectivity dis- 
tribution P r d{m) follows a Poisson law [2,13]: 



P rd (m) 



(3) 



For random networks, the connective constant can be 
obtained in a straightforward manner, due to the absence 
of correlations between links. For n = 1 one has obvi- 
ously ul d = z. Now, given a generic node and a link 
going out from it, the connectivity distribution Q r d(jn) 
for the other end of the link in a random graph is given 

by 



Then, the average number of two-step SAW's is given by 
u r 2 d = zJ2( m - l )Qrd{m) , (5) 

m>l 

and we find u r 2 — z 1 . Using the same procedure for 
n > 2, one has: 



(6) 



and thus for (large) random networks one finds /j, = z. In 
connection with this, we note that for a Bethe lattice (or 
Cayley tree) with connectivity z, the number of SAW's 
is given by u^ L — z(z — l) n_1 , and one has fi = z — 1 
(see, e.g. Ref. [38]). 



III. NUMERICAL SIMULATIONS 

The networks studied here have been generated from 
three different regular lattices: one-dimensional (ID) ring 
with coordination number z = 4 and 6, and 2D square 
lattice (z = 4). To construct our small- world networks, 
we consider in turn each of the bonds in the starting 
lattice and substitute it with a given probability p by a 
new bond. This means that one end of the bond (chosen 
at random) is changed to a new node taken randomly 
from the entire network. We impose three conditions: 1) 
no two nodes can have more than one bond connecting 
them, 2) no node can be connected by a link to itself, 
and 3) isolated sites (with zero links) arc not allowed. 
This method keeps constant the total number of links 
in the rewired networks (and consequently the average 
connectivity z). The total number of rewired links is 
\zpN on average. 

For networks generated in the present way there is a 
crossover size N* ~ p^ 1 that separates the large- and 
small- world regimes [39,40], and the small- world behavior 
appears for any finite value of p (0 < p < 1) as soon 
as the network is large enough. Our networks included 
1 x 10 5 sites in ID and 300 x 300 sites for the 2D system, 
so that we were in the small-world regime (system size 
N > N*). In the sequel we will call L the side length of 
the considered lattices, i.e. L — N 1 ^. Periodic boundary 
conditions were assumed. We note that our networks 
differ from those discussed by Watts and Strogatz [7] in 
that these authors left untouched z/2 links per site. For 
our simulated small- world networks we have obtained the 
average number u n of SAW's up to n = 21. Since u n 
increases with p, this maximum n was reduced to n = 14 
close to p = 1, in order to carry out averages over nodes 
of the generated networks. These numbers of steps in the 
SAW's are sufficient to obtain the connective constant ti 
with enough accuracy for our present purposes. In fact, 
the larger is p, the faster the ratio u„/u„_i converges 
with n. 
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FIG. 1. Mean-square end-to-end distance R„ for SAW's on 
small- world networks generated from a ID regular lattice with 
N = 10 5 sites and connectivity 2 = 4. We present values for 
i? 2 , normalized by the squared system length L 2 , as a func- 
tion of the number of steps n. Different symbols correspond 
to several values of the rewiring probability p. From top to 
bottom: p = 0.2, 0.1, 0.05, 0.02, and 0.01. A dotted line 
indicates the value R^/L 2 = 1/12 corresponding to p — 1. 
Dashed lines are guides to the eye. 



First of all, we calculate the mean-squared end-to-end 
distance of the walks on our small- world networks: 



i^ = <(r„-r ) 2 



(7) 



where ( ) indicates an average over n-steps SAW's with 
different starting sites and for different network realiza- 
tions with given p. Here, r„ refers to the position of site 
n in the d-dimensional Euclidean space of the underlying 
regular lattice, and r is the position of the origin for the 
considered walk. For SAW's on the ID lattices consid- 
ered here, one has R„ ~ bn 2 , whereas for the 2D square 
lattice ~ br?l 2 [29,35], with a lattice-dependent con- 
stant b of order unity in all cases. 

In Fig. 1 we present the ratio R„/L 2 as a function 
of the number of steps n for SAW's on ID small-world 
networks with z = 4, for several values of the rewiring 
probability p. R\ increases with n much faster than for 
the corresponding regular lattice due to the random con- 
nections introduced by the rewiring of links. (Note that 
for a regular lattice with L = 10 5 , we have for n = 20 a 
ratio R-n/L 2 ~ 10~ 7 .) In fact, if we call f n the fraction 
of n-steps paths that include at least one rewired link, 
we have: 



Rl « 6(1 - f n )n 2 + ^.fndL 2 



(8) 



The first and second terms on the r.h.s. come from SAW's 
without and with rewired links, respectively. The sec- 
ond term amounts, apart from the fraction /„, to the 
average Euclidean distance between any pair of sites in 



FIG. 2. Average number of self-avoiding walks u„ on 
small-world networks generated from ID rings with connec- 
tivity z = 4. We plot function of the path length n 
for several rewiring probabilities p, as derived from numerical 
simulations. From top to bottom: p = 1, 0.3, 0.1, 0.03, and 
0. Dotted lines are guides to the eye. 



a d-dimensional box with side length L. Thus, in the 
limit of large networks (as those considered here, with 
L ^> n) one has R^ w f n dL 2 / 12. In the course of our 
numerical simulations we have checked Eq. (8) by calcu- 
lating indcpendcnlty f n and R n as a function of n. We 
found that both sides of this equation coincide within er- 
ror bars (which are smaller than the symbol size in Fig. 
1). Then, R 2 n can be considered as a measure of the frac- 
tion of SAW's containing rewired links, and converges to 
dL 2 /12 for large n in the rewired networks (/„ — ► 1). 

We now turn to the number u n of SAW's on these 
networks. In Fig. 2 we show function of the 

walk length for networks built up from a ID lattice with 
coordination number z = 4. We have plotted results 
for several rewiring probabilities p, from p — (regular 
lattice, squares) to p — 1 (black diamonds). As expected, 
u n increases as p is raised, since introducing long-range 
connections in the starting lattice opens new ways for 
the SAW's. In particular, such long-range links allow the 
walks to visit regions far away from the origin for small n, 
and thus avoid the constriction associated to move close 
to the starting site, which limits the number of possible 
self-avoiding walks. In the logarithmic plot of Fig. 2 
one sees that log(u„/u„_i) converges rather fast to a 
constant for each rewiring probability p, which allows us 
to calculate the corresponding connective constant. Also, 
from the results shown in Fig. 2 we find that for large 
n the ratio between the number of SAW's for p = 1 and 
p = increases as k n with a constant k = 1.67. This 
means that at long distances, for each link available for 
SAW's in the regular lattice, we have in average 1.67 
connections for p = 1. This number is in fact the ratio 
between connective constants for p = 1 and p = in the 
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FIG. 3. Connective constant /iasa function of the rewiring 
probability p. Different symbols represent results obtained for 
small- world networks generated from ID regular lattices with 
connectivity z = 4 (squares) and z = 6 (circles), as well as 
from a 2D square lattice (diamonds). Dotted lines are guides 
to the eye. Error bars are less than the symbol size. Results 
for /i obtained by means of the analytical method described 
in Sec. IV are plotted as dashed lines. 



case d = 1, z = 4. 

The connective constant \x has been obtained for our 
simulated networks by finding the large-n limit of the ra- 
tio u n /u n -i. In Fig. 3 we present the resulting /i as a 
function of the rewiring probability p for our networks 
generated from ID and 2D lattices. One observes that 
\i changes fast close to p = 0, and the derivative dfi/dp 
decreases as p is raised. The largest change of \i in the 
whole region between p = and 1 is found for the net- 
works with z = 6. In fact, we find in this case an increase 
in fi of 2.13 to be compared with 1.48 and 1.07 in ID and 
2D rewired networks with z — 4. 

For p = 1 our numerical procedure gives in all cases 
connective constants p clearly lower than the mean con- 
nectivity z. In fact, we found /i = 3.69, 5.73, and 3.70 
(±0.01) for networks rewired from ID lattices with z = 
4 and 6, and from a 2D lattice (z = 4), respectively. The 
obtained values for networks with z = 4 coincide within 
error bars, irrespective of the starting ID or 2D lattice, 
indicating that for p = 1 the resulting rewired networks 
lost memory of the starting regular lattice. However, 
the \i values obtained for simulated networks with p = 1 
contrast with those expected for random networks, that 
coincide in each case with the average connectivity z, as 
explained above. This occurs because our networks with 
p = 1 are not true (poissonian) random networks, since 
they still keep memory of the starting regular lattices 
[19]. This memory effect is due mainly to the fact that 
one rewires only one end of each link, maintaining the 
other end on its original site. Hence, the connectivity 
distribution found for our rewired networks with p = 1 



P 

FIG. 4. Dependence upon the rewiring probability p of the 
change in connective constant, Afi — // — //o, with respect 
to the corresponding regular lattices. Symbols represent the 
same networks as in Fig. 3. Dashed lines were obtained from 
analytical calculations by using Eq. (13). 

does not coincide with that given above in Eq. (3). Such 
a difference with poissonian random networks should be 
even stronger for small-world networks generated in a 
way similar to ours, but leaving untouched z/2 links per 
site, as those studied in earlier works [7,21]. 

To analyze the change of \x as a function of p for a given 
underlying lattice, we call A/i = [i — /xo> Mo being the 
connective constant of the corresponding regular lattice. 
In Fig. 4 we show the obtained dependence of A\i upon 
p for the considered ID and 2D networks in a log-log 
plot. In all three cases we find that A^, can be fitted 
well by a power law A\i ~ p c for p < 0.01. For ID 
networks, the exponent c is found to be 0.98 ± 0.03 (for 
z = 4) and 0.99 ± 0.03 (for z = 6). For 2D networks we 
found c = 0.99 ± 0.03. Thus, our results indicate a linear 
dependence fj, — + a P for small p, irrespective of the 
underlying lattice. 

Therefore, the functional form for the p-dependence 
of A/j, close to p = docs not depend on the dimen- 
sion of the starting regular lattice. This constrasts with 
other properties of small- world networks, which have 
been found to change as p 1 ^, d being the dimension of 
the underlying lattice. This happens for example for the 
average number of nodes in "shell" n, which scales as p x l d 
[19]. Such a dependence is related to the scaling with d 
of the characteristic length scale of small- world networks, 
namely the average distance between ends of shortcuts £, 
given by £ = (pz)- 1 ^ [20]. 



IV. ANALYTICAL APPROXIMATION 

We now derive an approximate analytical expression 
that allows us to calculate u n in a small-world network, 
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assuming the sequence {u^} for the underlying regular 
lattice to be known. For a given path length n we obtain 
the mean number of SAW's u n by considering all possible 
sequences of unrewired and rewired links in small-world 
networks. With this purpose, we calculate the probabil- 
ity of reaching a rewired link as a function of the walk 
length. In particular, we will obtain the conditional prob- 
ability that the i'th link in a SAW is a rewired one, as- 
suming that link i — 1 is an unrewired one (i > 1). 

We first note that in a regular lattice the ratio Cj = 
Ui/u^_ 1 obviously depends on i. This ratio Cj measures 
the average number of available links starting from the 
(i — l)-th site in a generic SAW on the underlying regular 
lattice, and allows us to calculate the number of possible 
unrewired links in a SAW on a rewired network. Taking 
into account that the fraction of unrewired connections in 
the whole rewired network is 1 —p, then the mean number 
Qi of available unrewired links going out from site i — 1 
in a SAW is = Cj(l —p). On the other side, the mean 
number of rewired links going out from an arbitrary node 
reached by an unrewired link is given by (see Appendix 
A) 



Wi=p(z--). 



(9) 



(Contrary to Qi, the mean number Wi is independent of 
i.) Therefore, the conditional probability that link i(> 1) 
is a rewired one, assuming that link i — 1 is an unrewired 
one, is given by 



Pi 



Wi 



z p 



+ qi z'p + Ci(l -p) ' 



(10) 



with z' = z — 1/2. For i = 1, we take p\ = p. This 
probability pi is shown in Fig. 5 for the three types 
of networks considered here, for a rewiring probability 
p = 0A. 

Thus, the average number of i-steps SAW's that do not 
include rewired connections is 



(l-pi)...(l-Pi)«i. 



(11) 



and the number of those consisting of i — 1 unrewired 
links and a rewired one in step i is 



Bi = (1 - pi) ... (1 - pi-i) pi u° . 



(12) 



Now we note that a rewired connection in step i in a 
(large) small- world network ends on a random site of the 
network (most probably far away from the sites already 
visited in the same walk). This means that in step i + 1 
one begins most probably with a situation similar to that 
found in step i = 1. Finally, the average number of n- 
steps SAW's is given by 



^ ^ Bi 1 Bi 2 ...Bi j l Ai j 

i 1 + ...+ij=n 



(13) 



where the sum is extended to all possible combinations 
of indexes i\, with sum equal to n, including null 
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FIG. 5. Conditional probability p n that the n'th step in 
a SAW be a rewired link, assuming that link n — 1 is an 
unrewired one, as derived for p = 0.1 from the analytical 
procedure described in the text. p n is presented as a function 
of the step n for small- world networks built up from ID rings 
with z — 4 (squares) and z = 6 (circles) , as well as from a 2D 
square lattice (diamonds). Dotted lines are guides to the eye. 



indexes, for which we have Bo = Aq = 1. Each term in 

Eq. (13) represents a sequence of unrewired and rewired 

links, and thus the sum includes 2 n terms. The SAW's 

corresponding to the general term in the sum include j — 1 

rewired links (in steps number +t2,ii + ... 

As an example, Eq. (13) gives for n = 2: u 2 = B 2 + Bf + 

B 1 A 1 +A 2 . 

Note that the above equations, although rather accu- 
rate, are not exact. There are two reasons for this: First, 
each ratio Cj is an average number of allowed links start- 
ing from a site reached in i — 1 steps, but the actual 
number of such allowed links depends on the particu- 
lar site under consideration. Second, finite-size effects 
should appear, unless the considered networks are large 
enough. 

Values of connective constants fi derived from u n ob- 
tained with this procedure are presented in Figs. 3 and 4 
as a function of the rewiring probability p (dashed lines). 
There appears to be good agreement with \i derived from 
numerical simulations (symbols) for p < 0.2. For larger 
p the connective constants deduced from the analytical 
method are larger than those yielded by the simulations. 
(For z — 6 there is a region around p = 0.3 where the 
analytical results are slightly lower than those found for 
the simulated networks.) 

For p = 1 our analytical procedure gives \i = z, as 
for random networks. On the other side, our numerical 
simulations for small- world networks gave in all cases in 
the limit p = 1 connective constants /i clearly lower than 
the mean connectivity z, as indicated above and shown in 
Fig. 3. This difference between both procedures is due to 
the above-mentioned fact that simulated networks with 
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p = 1 do not have a poissonian connectivity distribution, 
which is implicitly assumed in the analytical method for 
this value of p. As a matter of fact, in this case we have 
Ai = for i > 0, Bi = for i > 1, and B\ = z, giving in 
Eq. (13) u n = z n , as for random networks [see Eq. (6)]. 
In this sense, our analytical procedure to calculate the 
number of SAW's gives a better interpolation between 
regular lattices and random graphs. 

In the context of this analytical approach, it is natu- 
ral to expect close to p = a linear dependence of the 
connective constant /j, on p, as found from our numerical 
simulations (see Sec. Ill), irrespective of the underlying 
lattice. By expanding the probability pi [Eq. (10)] to 
first order in p one finds for i > 1: pi = z'p/ci + 0(p 2 ). 
Introducing this expression for pi into Eqs. (11), (12), 
and (13), and keeping terms up to first order in p, we 
find: 

n 

u n = Y, B i A n-i (p<l), (14) 

i=0 

with the corresponding linearized expressions for Bi and 
A n -i. This expression includes contributions of SAW's 
containing zero (term A n , for i = 0) and one rewired 
links (all other terms, i = l,...,n). In this way, close 
to p = we find for the derivative d^i/dp the values 7.6 
(d = l,z = 4), 14.3 (d = l,z = 6), and 6.3 (d = 2). 
Thus, it is clear that the functional dependence of [i upon 
p for p <C 1 does not change with the dimension of the 
underlying lattice. 

V. CONCLUSIONS 

Self-avoiding walks provide us with an adequate tool to 
study the long-range characteristics of small-world net- 
works. For large networks, the number of SAW's in- 
creases asymptotically as u n ~ fj, n , provided that one 
considers system sizes L ^> n. For small-world networks 
generated from a given lattice, the effective connectivity 
\i ranges from the value of the regular lattice to \i = z 
for random graphs. For small p this effective connectivity 
follows a linear relation jj, = /xo + a P, a being a constant 
dependent on the underlying lattice. 

We have developed an analytical procedure to obtain 
the number of SAW's in small-world networks. This 
method is based on calculating probabilities of finding 
rewired or unrewired links in the walks, and gives re- 
sults in good agreement with numerical simulations for 
p < 0.2. The results of both methods differ for larger p, 
since they assume in practice different connectivity dis- 
tributions. Our analytical method gives in this respect 
a correct interpolation between regular lattices and ran- 
dom graphs. On the contrary, the rewiring (simulated) 
process gives rise to non-poissonian connectivity distri- 
butions, even for a rewiring probability p = 1, yielding 
in this case networks with connective constants /i lower 
than the average connectivity z. 



Both analytical calculations and simulations similar to 
those presented here can be useful to characterize other 
kinds of networks of current interest, such as scale- free 
networks, whose properties are known to depend on the 
asymptotic form of the connectivity distribution for large 
connectivities. 
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APPENDIX A: DERIVATION OF CONDITIONAL 
PROBABILITIES 

Here we calculate a conditional probability related 
to the number of possible connections starting from a 
generic node in a SAW, and necessary to derive Eq. (9) 
in our analytical approach in Sec. IV. In particular, we 
will obtain the average number of possible rewired links 
going out from a generic node X (so called for definitc- 
ness), assuming that this node was reached in a SAW 
through an unrewired link in step n. Here we mean by 
possible links all those connections that are available for 
step n + 1 in a SAW (leading to nodes not previously 
visited in the same walk). 

In a rewired network, the links with one end on node X 
can be classified for our present purpose into three types: 

1) Links which were not rewired and remain as in the 
original lattice. The probability distribution for the num- 
ber r of these connections is given by: 

Pi(r)=fy(l-p) r p z - r , r = 0,...,z. (Al) 

2) Rewired links for which the reference node X was 
not changed. Following the above notation, there are 
z — r rewired links, from which s keep one end on site X. 
The probability distribution for s is 



P2(S) = 




-r. (A2) 



3) New (rewired) links arriving at site X. The distri- 
bution of the number v of these connections is (for large 
system size N) 

P 3 (v) = ±- (tzy e- zt , v>0, (A3) 

with t = p/2. 

Thus, given a site with r unrewired connections, the 
number of rewired links is x = s + v , where s depends on 
r and v is independent of r. The probability distribution 
for x is: 

Q£\x)=Y^P 2 (s)P 3 (x-s), (A4) 

s=0 
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with s max = min{z — r, x). Then, we have 

*'«-ar--s('r)^^- 

(A5) 

Hence, the probability distribution for the number of 
outgoing rewired links (assuming that the incoming link 
was an unrewired one) is 

z 

Qout{x) = ]rQ m MQ^(x) , (A6) 

r=l 

where Qi n {r) is the probability of reaching a node having 
r unrewired links: 

Qin(r) = - z fy(l-p) r - 1 P z - r ■ (A7) 

[Qin(r) is given, apart from a normalization constant, by 
the product rP\ (r) .] 

Finally, the mean number of outgoing rewired links 
[calculated with the probability distribution Q ou t(x)] is 

oo 

(x) =^2xQ out (x). (A8) 

x=l 

Introducing expression (A6) into Eq. (A8), and after a 
straightforward but somewhat lengthy algebra, one finds 
the mean value: 

{x)=p{z-\). (A9) 

This average value (x) is called Wi in the main text, and 
the last result is Eq. (9) there. 
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